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ABSTRACT 

Stress intensity factor (K ) in a double cantilever beam 

ic 

with thin cantilevers is determined through a combined scheme of 
measuring strain near the crack tip and analysing data by using 
an FEM programme- The DCB specimen is made from two thin steel 
cantilevers bonded together with epoxy- The precrack is 

introduced during bonding by inserting a thin film of BOPP- 

The strain field of DCB specimen with thin cantilever cannot 
be expressed through a closed form solution because the free 
surfaces are close to the crack plane. Therefore, a finite 
element software has been developed for finding strain field 
around the crack tip for a stress intensity factor. This 

facilitates determination of K if strain close to crack tip is 

IC 

measured experimentally by establishing a relationship between 
strain and stress intensity factor. This relationship " is 

worked out by running the finite element code whenever the crack 
length is change- The finite element analysis also helps in 
selecting proper location for strain measurement. 

The specimen is loaded in mode I quasi -^stati cal 1 y on an 
INSTRDN machine under displacement controlled conditions. The 

strain is measured through a very small strain gauge of 0-2 mm 
gauge length, bonded at an angle of 45^ to the crack plane. As 

the crack grows and approaches the strain gauge, the output from 
the strain gauge increases till it reaches a peak- The peak 

corresponds to the location of the crack at about 1-5 mm prior to 
the strain gauge, measured along the- crack' plane. Further crack 
extension reduces the strain. The measured peak strain gives the 
value of critical stress intensity factor through the already 
established relationship between strain and SIF- The 

experimental 1 y obtained value of has been compared with 

numerical values obtained through finite element analysis. The 
measured value of K is found to be lower. The reason is 

IC 

thought to be imperfect bonding between the cantilevers due to 
air entrapment. 
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CHAPTER I 
INTRODUCTION 


In composite materials, one of the major concerns is the 

poor interlaminar toughness on impact. Under impact loading the 

interlaminar cracks can propagate at high velocities exceeding 

300 m/s, according to observations made using high speed cameras 

C13. The problem to be investigated is the relationship between 

fracture toughness and the crack velocity. This information 

plays an important role in improving the performance of composite 

materials. The toughness of interlaminar crack growth may be 

measured in terms of Critical Energy Release rate (6 and G ) 

ic lie 

or Critical Stress Intensity Factor (K and K >. 

IC lie 

1.1 MEASUREMENT OF 6 AND G 

IC lie 

The method of energy balance applied in the measurement of 

dynamic 6^ or 6^^ is inaccurate as it is not always possible to 

account for all the energy losses since dynamic experiments deal 

with wave propagation inside the specimen. As wave propagates, 

it interacts with defects causing transformation of kinetic and 

strain energy into heat. Such transf ormation of energy is also 

known as deformation friction C23. The heat produced may 

increase the temperature of the specimen or part of it may be 

transferred to the specimen holder or surrounding air. As this 

heat generation cannot be measured accurately, it affects the 

results substantially since G for interlaminar fracture in 

2 

composites is very small, of the order of 200 J/m , compared to 

2 

metals where G is more than 30,000 J/m . Hence, we find that 
experimental determination of dynamic G or G is difficult and 

1 C I IC 

the stress approach is more suitable. 
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1.2 MEASUREMENT OF K OR K 

ZC lie 

For specimens with thickness of the order of uncracked 
ligament, equations are available for calculating stress/strain 
field in the specimen through the stress intensity factor. But 
for very thin specimens no such closed form solutions are 
available in the literature to the best of author's knowledge. 
Since in the current analysis we are using thin laminate, free 
lateral surfaces are close to the crack tip (<3 mm thick) and 
stress or strain fields cannot be predicted by conventional 
formulae. Moreover, wave propagation effects become prominent 
under dynamic conditions. So the measurement of K is not 
straightforward here. 

The strategy adopted is to solve the problem in stages. In 
the first stage, a double cantilever beam (DCB) specimen made of 
thin isotropic metal (steel) strips bonded with epoxy will be 
tested for fracture toughness. Assisted by finite element 
analysis, a proper location will be selected for fixing a small 
strain gauge to measure the strain near the crack tip in mode I 
under quasi— static conditions. Based on the results of the 
finite element analysis, a relationship also will be established 
between strain and stress intensity factor. A strain gauge of 
0.2 mm gauge length will be bonded at the location and the DCB 
specimen will be loaded in mode I. The measured strain will give 
the value of through the relationship established between the 
two. In the second stage, following the same line, use of strain 
gauges on isotropic materials under dynamic conditions will be 
investigated. Finally in the third stage, with enough confidence 
derived from the above two stages, strain gauges will be used on 
fibre composite materials in dynamic experiments. In this study, 
work has been done only on the first stage. 
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1.3 LITERATURE SURVEY 


Work on stress/strai n field in a double cantilever beam with 
thin cantilevers loaded statically has not been so far reported 
in the literature to the best of author's knowledge- However, 

expression for K of such a specimen is listed in E83. This is 

Id 

derived in an indirect way by first finding compliance, C, of the 


specimen and then finding G by the formula G = ~ ^ 


2 B 

P is the load on the cantilever, a is the crack length and B 
the width of the specimen- 


where 
i s 


is then determined by” using the 


formal i 


G E / (l-h> y 

X 


The proof is given in Appendix B- 
In the derivation only flexural energy of the cantilever has been 
accounted for; that is , the shear strain energy caused by the 
lateral loads at the cantilevers and the strain energy in the 
bulk material beyond the crack tip have not been considered in 
evaluating the compliance- However , the assumption of neglecting 
these strain energies is reasonable because the cantilevers are 
thin and flexural energy dominates. Another point to be noted is 
that the procedure only provides the value of but gives no 
idea of strain field around the crack tip which is strongly 
influenced by the nearby traction free surfaces (Sec 2.8-1)- 


Since work on stress/strai n field in DC5 specimen is net 
available in literature, this survey on literatuf“’e is presereted 
on two related fields viz-, (i) stress field around a crack tip 
in a large FRP plate and (ii) experimental evaluation of dynamic 
in large metal or plastic plates- 


Rosakis, Duffy and Freund IZl performed dyna.mic crack 
propagation experiments using wedge loaded double cantilever 
beam specimens of an austeni tized , quenched and tempered 43^D 
steel- They measured the dynamic stress intensitv factor by 
means of the optical method of caustics. The e<peri meet a.l data 
was obtained from the shadow spot patterns photographed with a 
Cranz— Schar'di n high speed camera- The i nterpretati on of the data 
was based on an el astodynami c analysis. The- obtained 


the 
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instantaneous value of the dynamic stress intensity factor as 
a function of crack tip velocity. It is observed that the 
fracture toughness increases with increase in crack velocity. The 
specimen used by them was large in size and therefore the 
technique cannot be used for DCB specimen with thin cantilevers. 

Stress Intensity Factor of a crack in a large orthotropic 
composite materials was measured using strain gauges by Bhukia, 
Agarwal and Bharat Bhushan C4I. They first developed theoretical 
equations for the strain field in the vicinity of the crack tip 
in an orthotropic material. These equations were then critically 
evaluated to obtain the optimum location and orientation of the 
strain gauges to be used for strain measurements. They then 
verified the appi i cabi 1 i ty of the procedure using a glass epoxy 
specimen. This work deals with large plate only and that too 
under quasi— static loading conditions. 

Strain gauges were used by Berger and Daily to measure the 
dynamic fracture toughness of steel 153. They recorded a 
sequence of strain time traces during the passage of crack 
beneath a line of strain gauges. The data analysis uses the 
spatial distribution of strains along the gauge line at a 
specified time. They developed an overdetermined system of 
equations which are solved to determine K as a function of time 

ID 

and position during crack propagation. In this study also, only 
large plates are dealt with. 

r4ishiDka et. al « LSI studied the specimen size effects oi! 
dynamic crack propagation and arrest in DCB specimens using laser 
caustic method. They worked on polymethyl methacrylate (PMHA) 
made specimens in mode I. Dynsttnic condition was obtained using 
the conventional technique of storing energy through a blunt 
crack. They worked with maximum crack velocities of the order of 
200 m/s« The specimen length was changed in order to study the 
specimen size effects on dynamic crack propagation and arrest- 
They found that dynamic stress intensity factor decreases as the 



1 ncreases 


crack moves forward. They also observed that K 

X0 

crack velocity- However the studies were conducted on 
speci mens. 


wi th 
1 arge 


Kolednik in his theoretical study £73 gave a physical 
interpretation of the J— R curves for elastic-plastic fracture- 
He derived the difference between the crack-i ni ti at i on toughness 
and crack-growth toughness using energy balance under 
quasi -static conditions. He considered three point bend 
specimens each of which consists of two parts glued together 
along the ligament. Here also large specimens were analyzed- 


1.4 WORK DONE IN THIS STUDY 


A DCB specimen has been prepared which simulates the 
delamination in composite materials, when loaded in mode I. It 
is made fay gluing two thin hardened alloy steel strips together 
with epoxy- A finite element code has been developed for 
estimating the strain field near the crack tip in the specimen 
which was then related with stress intensity factor- The strain 
estimations have been made use of in selecting a proper position 
and orientation for fixing the strain gauge- The specimen has 
been loaded in mode I (Fig- 1-1) under displacement controlled 
conditions. The strains near the crack tip have been measured - 
This observation is matched with strain from the finite element 
programme and the value of corresponding to the observed strain 
is obtained- 


The finite element formulation is dealt with in chapter II 
and the experimental technique is detailed in chapter III. The 
results are discussed in the IV chapter- 
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P 


t 



Fig, 1.1 A DCB Specimen. 
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CHAPTER II 

STRAIN FIELD NEAR CRACK TIP OF A DCB 
SPECIMEN IN MODE I THROUGH FEM 


2.1 INTRODUCTION 

In this chapter, numerical analysis has been carried out 
through FEM to determine the strain field near the crack tip for 
a double cantilever beam specimen with thin cantilevers, in terms 
of loads applied at the cantilever ends- The goal of this 
simulation is to identify a suitable location for a strain gauge 
near the crack tip so as to optimize its output and relate it to 
stress intensity factor. 

The geometry of the DCB specimen used is shown in Fig. 2.1. 
It is made by joining two thin strips of hardened alloy steel 
with epoxy. The epoxy used is LY 556 and the hardener is HY 1907 
IN. The curing is done at a temperature of 130 degree Celsius and 
under a pressure of nearly 4.9 MPa for one hour. The thickness 
of the epoxy bond is very small, of the order of 0.03mm and is 
considered to be a line in the numerical formulation. However it 
is to be noted that the fracture behaviour of the specimen will 
mostly be determined by the bond. 

2.2 SPECIMEN BEDMETRY AND MATERIAL 

* 

The alloy steel used is hardened 40Ni2CrlMQ2S <Er4 24). 
Hardening increases its yield stress value thereby dec?^easing the 
size of the plastic zone near the crack tip, which is estimated 
by Irwin's formula C81, 


r 

p 


ys 



2.78mmT 



CRACK LENGTH 
VARIED BETWEEN 
25 and 50 mm 


Fig. 2.1 Test Specimen. 
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In the present case, specimen has a Rockwell hardness o-f 33 
on the C-5cale and a yield stress of around 1000 MPa- In order 
to estimate the size of the plastic zone, of the DCB specimen 
was determined through some preliminary experiments by finding 
G and then using the relation K - JG E~/ for plane 
strain. For of the order of 1.5 MPa-^m, the plastic zone size 
3S less than ©.0007mm which is small enough to be neglected and 
the whole region can be considered to be elastic. 

2.3 FORMULATION 


In this 2-“D finite element formulation, the standard method 
of displacement based stiffness formulation is used C93- The 
technique is widely used but for the sake of completeness, it is 
briefly explained below- Since the dimension along the third 
direction (Z— direction of Fig- 2.1) is large and since the hinges 
used to apply load distributes the load uniformly in the 
2~direction (Fig. 1.1), this is a plane strain problem. 


The objective is to discretize the domain into a number of 
finite elements, assume a displacement variation in each, 
construct the potential energy § and minimize § of the whole 
body. The potential energy is defined as 


§ = U - N 

ext 


. . . ( 2 . 1 ) 


where U represents the strain energy and » the; work done 
by external forces. External forces include boundary forces and 
body forces. In the present analysis, body forces have been 
negl ected . 


For a linear elastic material, the strain energy pe? unit 
volume is given as 

dU = r cy d£ = - cy . . . Cr.Z) 

j L) tj 2 tj I.J 


O 
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where are stress and strain components at a point 

within the element- 


Thus, the strain energy of the element in a plane problem 
becomes , 


^ - r a £ - rr <y £: o r ^ 1 

o 2 ^ ij ij 2 J *- XX XX yy yy xy xy -* 


dV 


Using the matri)? notations for stress and strain components, 


= < <y 


yy 


xy 


< 2 . 3 ) 


= < s 


yy 


xy 


<2. 4) 


The strain energy can be expressed in the fo^-m, 

£§•3 = - r i£y^€<y} dV ...(2.5? 

e 2 J 

V 

where v represents the volume of an element. 

Uork done by surface forces is 

r r F u dS ... (2-6) 

€'Xt J 

s 

where u is the displacement, F is the force and s is the surface 
on which F acts. For the 8~noded i soparametr i c element chosen 
(Fig. 2.2), the di spl acements are given by 


X 

u = a + 
o 

a X 

1 

-t~ a 

2 


2 

a X 

3 

a X y 

a 

5 

2 

V 

4- o 

<5 

2 ^ : 
X V 4- a -iy 

f 7 

= b + 
o 

b ): 

1 

b 

2 

y -*■ 

b x^ 

3 

+ b X y 

4 

b 

5 

z 

y 

4- b 

<5 

2 I 

X y 4- b xy 

7 


. . (2.7) 
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where 


u 


1 

X 

i 

2 


u 


X 

U , U 


at point (x ,y ) 
11 

at point <x ,y ) 
2 2 


X y 

u , at point (x ,y ) 

0 6 O fi 


Alternatively, u , u can be expressed in terms o-f 


nodal 


displacements in the form 


X Iwl ILI ^ . 

u = Nu+Nu + 
11 2 2 

= N u^+ N u^+ 


N u 
6 8 


E N 


t L 

y _ f-. ki ..y 


1 1 


2 2 


-I- N = r N u 

6 8 ^ i L 


1 


1 


. . . <2.B> 


The geometry mapping can be similarly given as 


X = r N X 
t t 

y = E N y 

“ V V 




. . ( 2 . 9 ) 


where N ,N etc. are functions of x,y,x ,y ,x ,y . . . . 
12 112 2^ 

Putting the above two equations into matrix form. 


'V^e- 


r X -1 

U 

y = ^ 


0 

N 

z 

0 . . 

N 

a 

0 ■ 



[ 0 

N 

1 

0 

N 

2 

0 

N 

B-* 


u 

u 

u 

u' 

2 

y 

U 

G 

\ 

U 


8 


= CN3 IlO 


. . . (2. 10) 
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The strains can be obtained by differentiating the displacements 
as 


s 

XX 


■ au^ 

ax 

« jT* 


du^ 

yy 


dy 



gu%guy 

y 

xyj 


ay ax 


which can be further expressed as. 


m m 


dN 


aN 



s 


1 

0 

2 

C5 

. 0 

XX 


ax 

dx 






dN 


dN 

dN 



0 

1 

0 

2 

8 

yy 


ay 

ay 

dy 



dN 

aN 

dN 

aw 

dN 



1 

X 

2 

2 

8 

^xy J 


ay 

ax 

dy 

dx ------ 

dx 


u 


u 


u 


I- B -» 




... ( 2 . 11 ) 


The stress strain relation is represented as 

io>y = CDl Cc> 

where CD! is the material property matrix defined as 

01 


CD3 


(1+1^) (l-2i>) 


1-v V 
V 1—v 

0 0 


0 

l-2v 


for plane strain. 


In the above matrix, v is the Poisson's ratio and E is the 
Young's modulus. Now, the expression for strain energy can be 
rewritten in the form. 


U = I dV = - J{:uj'*’CB3'*‘CD3i:B3Cu>t dxdy ... (2.12) 


V 


A 
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where t is the thickness of the element. Transforming into 
natural coordinates f and t). 


I Jfu>'*’CB]’‘CDDEB3£u>t |J I d?dT) 

A 


where |J| is the Jacobian. Substituting in Eq.2.1, 


. . . (2. 13) 


t J fu>'*'CB3’'CD3EBD<:u>t |Jj d?d7> - J F u dS 


- fu>'*'cK®3fu> - Cu>’’CF®> 
2 


<2. 14) 


where 

CK3 = J CB]’'eD 3CB3 |J| d? dr) 

A 

and the superscript e represents an element. The total potential 

energy of the system can be obtained by adding up elemental 

potential energies. 

N N N 

§ = j; = £; CuT EK®3 Eu3 - E fuJ -CF^J 

essl e = l 




- 


J. a,T---a- 

•tu > -CF > 


. . . (2. 15) 


Where the superscript o denotes the global domain. It is 

N 

O & 

noted here that K = E represents assembly and not 

©= 1 


addition. 


to be 
direct 


To find the displacement u which minimizes §, is equated to 

zero. From equation 2.3, 

|t = CKJ-CU^J - (F°y = 0 
du 

i.e., EK3°Tu°> = £F°3 ...(2.16) 


The set of simultaneous equations EK3£u> = -CFJ has been 
solved using a Sky-line solution algorithm. 
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In the analysis, the load is applied at the end of the 
cantilevers so that opening mode conditions are simuiated« 

2.4 MESH BENERATION 

The mesh used is shown in Fig. 2-3 and 2.4. The former 
shows the fine and coarse elements used whereas the latter shows 
the collapsed elements at the crack tip. The elements are 
8~noded isoparametric ones- Near the crack tip, elements are 
kept fine and coarse elements are chosen away from the crack tip 
where stress variations are small. The size of the small 
elements which lie near the crack tip is 0.5 mm. The coarse 
elements have a size of 13 mm- 

At the crack tip collapsed quarter point elements are made 
use of to model the square root singularity. The behaviour of 
quarter point element is explained in Appendix A- In our 
analysis collapsed quarter point elements are used which has the 
said property along all directions towards the singular point 
C113. 


2.5. CALCULATION OF STRESS INTENSITY FACTOR 

The stress intensity facto^^ can be determined as follows 
191 (Fig. 2.5) 

K = £_£ I r _Il ') r (-4 V - V 1 - ("4 V - V ]1 ... (1. 17) 

I Q+i ^ 2L ^ B2 C2-' B1 cl-' J 

Where Q = 3 — 4i^ 

G = Sheer Modulus 


and 


u = x-dispi acement at the node point 
V = y-di spl acement at the node point 




Fig. 2.3 Mesh Showing Coarse and Fine Elements for 
DGB Specimen with Thin Cantilevers. 



Fig. 2.4 Collapsed Elements at Crack Tip for 
DCB Specimen with Thin Cantilevers 
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Fig. 2.5 Collapsed Element Near Crack Tip. 
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2.6 FLOW CHART 


The flow chart for the programme is as shown in Fig. 2.6. 



Fig. 2.6 Flow chart for finite element programme 

The software has been run on HP 9000 and CONVEX 220 systems. 
The latter is a mini super computer approximately six times faster 
than the super mini HP 9000. The execution time for a mesh of 
1384 nodes with 364 elements is nearly 30 seconds on CONVEX. 

2.7 TEST PROBLEMS 

The programme has been checked using following test problems 
of which results are known. 

1. Large plate with a center crack 

2. Large plate with an edge crack 

2.7.1. LARBE PLATE WITH A CENTER CRACK 

The problem of large plate with center crack under 
distributed load is solved for stresses and stress intensity 
factor in Mode I. Figure 2.7. shows the geometry of the specimen 
and Fig. 2.8‘, the mesh used. 
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The FEM code gave a stress intensity factor of 3.57 MPa-^m 

compared to a value of 3.76 MPa^m' f rom the formula K»=J Sec o-J tta 
C83. The two are in close agreement, the numerical result 
differing from the latter by 5 % only. 


Also in this case of large plate, the stress intensity 
factor gives the stresses near the crack tip through the 
equations given below CSl : 


a 

X 


a = 
y 


T 

xy 
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The stress field obtained numerically and found to be in good 
agreement with analytical results. 

2.7.2. LARGE PLATE WITH AN EDGE CRACK 

The geometry of the problem is shown in Fig. 2.9. The 
stress intensity factor and stresses in opening mode are 

solved. Figure 2.10. shows the mesh used. 

The value of K is obtained as 0.32 MPa-Jm from finite 

I _ 

element calculations. This is within 3 X of the value 0.33 MPa-|m 
calculated using the equation K = 1.12 o^ttaL . Here also, the 
stress field is found to be in good agreement with analytical 
results. 


2.8 RESULTS OF DCS SPECIMEN WITH THIN CANTILEVERS 

Finally, a DCB specimen of present investigation is analyzed 
for stresses and SIF. The specimen geometry is the same as in 
Fig. 2.1. The mesh is shown in Fig. 2.3 and 2.4. 



200mm 


BOON 



Fig. 2. 9 Problem of Edge Crack 
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The value of stress intensity factor is obtained as 
1.6 MPa-Jm from finite element formulation for a load of P = 60N. 

2.8.1 Comparison of Numerical Result 


•^o^P^i^ison purpose, has also been calculated using 
the formula C83 


K 

I 


2-/3 P a 

. 3/2 2 . 

h B ( 1 -p ) 


. . . (2. 18) 


This method is referred to as ' through approach ' in 
this study since Eq. 2.18 is derived in an indirect way, first 
finding G and then using the relation K = I — , The 

derivation of Eq. 2.18 is presented in Appendix B. 


This analysis provides by knowing load at the cantilever 
ends, crack length and specimen dimensions. However, it does 
not give any idea of strain field near the crack tip. In fact 
obtaining closed form solution for strain field is difficult 
task since the free surfaces are very close to the crack tip 
and influence the stress/strain field considerably ever. _L 
points which are close to crack tip. Because the strain gauge i« 
of finite dimensions (0.2 mm gauge length and 1.2 mm width), 
it cannot be fixed very close (< 1 mm) to the crack tip and 
for strain field in a DCB specimen with cantilevers of 
thickness less than 3 mm, one is forced to go to the finite 
element analysis. 


Moreover, it is clear from the derivation that only flexural 
energy of the cantilevers has been accounted for and the shear 
strain energy caused by the lateral loads and the strain energy 
in the bulk material beyond the crack tip have not been 
considered in evaluating the compliance, although, the assumption 
of neglecting these strain energies is reasonable since the 
cantilevers are thin and flexural energy dominates. 
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Using Eq. 2.18, the value of is obtained as 1.89 MPaym 
for load P = 60 N. This differs from the result of finite 
element analysis by IS X. The numerical programme has been run 
for different loads and crack lengths. The results are compared 
with those calculated through Eq. 2.18 in Table I. The 
percentage difference between the values obtained through the two 
methods is consistent. 

For comparing the experimental results the numerical value 
is preferred because it is free from assumptions made in the 
approach. 


TABLE I Comparison between K obtained numerically and K 
calculated through equation 2.18 


SI . 
No. 

Load 

N 

Crack 

1 ength 

mm 

through 

G approach 
^ MPaym 

through 

F.E. analysis 
MPaVm 

diff- 

erence 

7. 

1 

58 

36.5 

2.29 

1.89 

-18 

2 

50 

36.5 

1.98 

1.63 

-18 

3 

56 

36.5 

2.21 

1.82 

-18 

4 

70 

35.5 

2.77 

2.29 

-17 

5 

82 

30.5 

2.71 

2. 18 

-19 

6 

81 

23.5 

1.93 

1.76 

-14 


2.8.2 Strain Field and Contours 

The results are post processed to obtain the strain field 
around crack tip in the thin specimen of Fig. 2.1 using a routine 
developed by Ramesh E121 and strain contours are plotted. 


First step while plotting is to solve for the strains at 
Gauss Points in each element. These strains are then linearly 
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extrapolated to find strains at the nodes. For extrapolation, 
the following formula is used C133. 
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Where e , b , s ^ e and c, ^ s. , ^ e, . . . o, each represents one 

A B C D l ll lll lV 

strain component e , e or ^ . The roman numerical subscripts 

X y xy 

represent the nodal quantities and the alphabetical subscripts 
represent the gauss point quantities. In 2 point integration, p 
= (Fig. 2.11). The principal strains and their orientations 

are then determined at each node. For finding the principal 
strains, the following formula is used C143 


£ 

1,2 



The orientations of the principal strains are calculated by 
solving for eigen vectors of the strain matrix corresponding to 
£^ and £^ which are eigen values. 

2.8.3 Results 

The strain fields are given in Fig. 2.12 to 2.14. All the 
plots are for the upper cantilever. Fig. 2.12 gives the strains 
along the y-axis and Figure 2.13 shows the strains along x-axis. 
It is clear that the strain values are high near the crack tip. 
The principal strains are also plotted which is shown in 
Fig. 2.14. 
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2»8-4 Location of Strain Gauge 

It is noted that the points which lie on the line of 
intersection of the neutral plane of the upper cantilever and the 
side face of the specimen and close to the crack tip has 
relatively high strains suitable for measurement and does not 
- ’ 7 '^*^'* ence any bending strain after the crack tip passes the 
point- Of these points, the one which is at a distance of 
approximately 1-5 mm from the crack tip (measured parallel 
to the x-~axis) experience maximum strain 

At this point, the orientation of principal strain is 

approximately 4S^ with the x- axis (Fig- 2.15) and hence it is 
decided to fix the strain gauge at a convenient angle of 45*^. 
Figure 2.16 shows the strain contours along 45^ direction. The 

location of strain gauge and the area covered by it are shown in 
Fig- 2-17. 
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Fig. 2.12 ^yy for DCB Specimen with Thin Cantilevers 




Fig. 2.13 for DCB Specimen with Thin Cantilevers 
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Fig. 2.14 Principal Strains for DCB Specimen with 
Thin Cantilevers 




Fig. 2.15 Orientations of Principal Strains for DCB Specimen 
with Thin Cantilevers 
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Fig. 2.16 Strains along Direction that is at 45* 
with X-aixs. 
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Fig. 2.17 Strains along 45'’ with Strain Gauge Location, 


CHAPTER III 


EXPERIMENTAL TECHNIQUE 


3.1 INTRODUCTION 

This chapter deals with experimental details including 
design and preparation of specimen, equipment used, parameters 
controlled and measurements made during the experiment. The 

experiment is conducted under controlled displacement condition 
so as to obtain a quasi— static crack growth. 

3-2 SPECIMEN GEOMETRY AND MATERIAL 

GEOMETRY 

The DCS specimen used has two thin metal strips joined 
together with epoxy. The crack length is varied between 23.5 and 
36.5 mm and the thickness of the cantilevers is 2-78 mm each. 
Figure 2.1 provides information on the geometry. 

MATERIAL 

The material used is hardened alloy steel, 40Ni2CrlMo2S 
(EN24) , which has a high yield stress value after hardening. 
Hardening is done by keeping the metal at 830 — 850 '°C for one 

hour and then quenching in oil. This is followed by annealing to 
relieve the residual stresses. In annealing, the metal is heated 
to about 310*^0 for a period of nearly 20 minutes accompanied by 
furnace cooling. The hardness of the strips is then checked. In 
the present case the average hardness of the strips on C-scaie is 
33 which corresponds a yield stress of nearly 1000 MPa. 
strip was ground flat and polished with alumina powder. 


Each 



BONDING 


The "faces 0 "f the metal strips to be bonded using epoxy are 
etched to give better inter -facial strength. The etching was 
done by dilute Nitric Acid. For bonding, epoxy LY 556 and 
hardener HY 1907IN in the ratio 100:85 by weight are used. To 
this mixture, accelerator DY 062 also is added at the rate of 1.5 
ml per 100 gm o-f epoxy. The resin, the hardener and the 
accelerator are supplied by Ciba-Beigy Limited, Bombay. The 
sides ot the specimen are covered with a layer of wax to prevent 
them from getting coated with epoxy and the faces to be bonded 
are cleaned with pure acetone. The mixture is then applied to 
the cleaned faces. The two strips are then put together and 
pressed in specially built fixture made of two stiff steel plates 
(each 25mm thick) with ground flat faces (Fig. 3.1). Between the 
strips and the flats of the fixture, a sheet of BOPP was placed 
on each side that works as a release film. 

Provisions are made in the fixtures for keeping the two 
strips in position while curing. The base plate of the fixture 
has two cylindrical stops and three adjustable buttons. Both the 
thin strips are butted against the cylindrical stops and kept in 
place with the adjustable buttons. Each adjustable button has 
slot of 4 X 7 mm so that they can be tightened to the base at a 
desired location. Also, these buttons have different heights so 
as to hold each strip in its place. 

The fixture along with the specimen is placed between two 
platens of a hydraulic press. The platens are heated through 
inbuilt heaters and the temperature of the specimen is monitored 
by placing a chrome! -al umel thermocoupl e. The temperature of the 
specimen is maintained at 130*^ C under 1.2 MPa pressure for one 
hour. It is then allowed to cool under the same pressure to room 
temperature. 



PRECRACK 


While bonding the strips, a crack is introduced by placing a 
•film of Bi axially Oriented Poly-Propylene (BOPP) up to the 
required length so that the strips are not bonded there. For 
experiments of this study, crack length is varied between 23.5 
and 36.5 mm. 

3.3 HINGES 

Hinges have been designed and fabricated for applying load 
to the cantilever ends (Fig. 3.3). For fixing the hinges, two 
tapped holes have been provided at the ends of the cantilevers as 
shown in Fig. 2.1. 

Each hinge consists of a bracket, a support and a pin all 
made of steel. The bracket is fixed to the cantilever with the 
help of two high strength capscrews. The support is inserted 
into the bracket and the pin is used to join the two. The 10 mm 
diameter pin makes the hinge strong and stiff. All the holes of 
the hinges and the diameter of the pin are made with close 
tolerances. The end of the support is cylindrical of diameter 
12mm, made to fit the inside loading fixture of the tension 
machine, INSTRON 1195. 

The hinge so designed is moment free while load is applied. 
Further, it keeps the load line always vertical. 


3.4 EXPERIMENTAL DETERMINATION OF G 


ic 


A preliminary experimental study was conducted to have some 
idea about the maximum load as well as order of S through the 
standard technique of successive growth of the crack. The details 


of the technique are given in C153. Knowing 


1C 


IC 


IC 


one 


can 


G E/ ) - 

IG 


evaluate K 


|iy usi ngi the well knovvn relB-tion K 
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3.5 EXPERIMENTAL DETERMINATION OF K 

rc 

3.5.1 Introduction 

This constitutes the major objective in the present work. 
The idea here is to correlate strain measurements in the vicinity 
of the crack tip through a strain gauge with stress intensity 
factor for a DCB specimen with thin cantilevers. The loads are 
applied in displacement control mode with a very slow cross head 
speed. This work is expected to pave way for a more widespread 
use of strain gauges in this field. It is hoped that this 

procedure could be extended to dynamic problems also. 

3.5.2 Strain Gauges 

The strain gauges used in these experiments have been 
supplied by Tokyo Sokki Kenkyujo Co. ,Ltd. , Japan. They have a 
gauge length of 0.2 mm and a gauge width of 1-5 mm. The base of 

the strain gauge is rectangular with dimensions 3.5 mm X 2.5 mm. 

Its resistance is 120 +0.3 O and gauge factor, 2-07. Figure 3.4 
shows the details of the strain gauge. 

Fixing of Strain Gauges 

Extreme care is taken in fixing the strain gauges as it 
has only a very small area for bonding. Epoxy LY 556 and 

hardener HY 951 are used in bonding. The two are mixed well in 

the ratio 10:1 by weight. The ambient temperature is kept 

around 55 to 60‘^C for about 24 hours to help curing. Localized 
pressure is applied while the curing takes place. 

Positioning of the strain gauge is done with the help of a 
fixture made out of self sticking paper. It gives a corner and 
two sides as reference for fixing the strain gauge (Fig. .-^.5)- 
The efficiency of this method is checked from the observations 
made under travelling microscope which showed the center of the 
strain gauge to be very near to the point required, the error 
being less than 0.1mm. 



40 


RELEASE 



Fig. 3.1 Pressing of Two Steel Cantilevers for Bonding them with 
Epoxy at High Temperature and Pressure with a pre-crack 
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3.5.3 Measurement of Strain 

The strain gauge output is measured using strain indicator 
(Model P-350A, Measurement Group). The crack velocity being very 
small the strain indicator can be conveniently used. The strain 
indicator has a least count of 1 p and hence provides reasonably 
accurate observations. The measurement is made continually at 
regular intervals o-f time. 

3.5.4 Equipment 

The experiment is conducted on INSTRDN, model 1195. This is 
a highly flexible machine and can be used for testing materials 
in tension, compression, flexure and torsion. This has two main 
units namely an electronic control console and a loading frame. 
There is a chart recorder system also which is used in our 
experiments to plot the displacement against load. INSTRDN has a 
maximum load capacity of 100 kN and a lowest cross head speed of 
0.05 mm per minute. The cross head speed can be maintained 
constant with an accuracy of +0.1 7. when averaged over 20 mm of 
travel at constant load. 

3- 5- 5 Procedure 

The specimen is mounted on the INSTRDN machine and the 
trailing end of the specimen is suspended by ' t-'.. csaO to keep the 
axis of the specimen horizontal- It is a pulley like arrangement 
so that the specimen as a whole can move up and down but remains 
always horizontal. 

The specimen is loaded in di spl acement controlled mode and 
the cross head speed is fixed at 0.1 mm/minute. From the 
preliminary experiments it was observed that the crack velocity 
remains quite low under such a displacement rate. Since the 
maximum load during the initial experiments was in the range of 
50 - 75 N, the full scale load on the machine is kept 100 N. The 
chart speed for recording the graph displacement versus load is 
maintained at a convenient value of 10 mm/minute. In this 
machine, the upper jaw is fixed and the lower jaw moves down thus 
applying the load. 
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Once the experiment is started and the load reaches around 
halt the expected maximum loady recording of the strain 
indicator readings is started. Once the crack has moved well 
past the strain gauge, the specimen is unloaded and the readings 
are noted while unloading also. The entire experiment lasts -for 
approximately 20 minutes. 

3.5.6 Data analysis 

The experiments were conducted with the strain gauge fixed 
at a suitable distance from the center line of the hinges so that 
it is at about 5 to 8 mm from the tip of the precrack- The 
smaller value is used in the case of smaller crack lengths. From 
the numerical analysis, it is estimated that the maximum value of 
strain recorded will be when the crack tip is 1.5 mm away from 
the strain gauge. Hence, the crack length corresponding to the 
peak obtained in the strain plots can be deduced. We find that 
for a crack length of 30 mm and load of the order of 60 N, the 
maximum strain expected is around 140 microns according to 
numerical prediction. The graphs and other results of the 
experiments are shown in subsequent pages. 

3.6 CHECKING DF STRAIN GAUGE 

The performance of the strain gauge of 0.2 mm gauge length 
has been checked by measuring the bending strain in a cantilever 
in static condition. The error in measurement was found to be 
within 5 %. 




Fig. 3.2 Photograph of Bas 


Fixture used in Cunng 





Fig. 3.4 Detailed Figure of Strain Gauge. 
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Fig. 3.5 Fixture for Strain Gauge through Self Sticking Papers. 
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CHAPTER IV 

RESULTS AND DISCUSSIONS 


4.1 INTRODUCTION 

The strain field in the DCB specimen with thin cantilevers 
has been obtained through finite element technique and used in 
selecting proper location for fixing a strain gauge for strain 
measurement as described in section 2.8.3. After gluing the 
strain gauge of 0,2 mm gauge length in the ’proper position (Fig. 
2.18) , the specimen is loaded in a tension machine in Mode I 
at a very low crosshead speed. The load applied at the ends of 
the cantilevers is recorded through a strip chart recorder of the 
tensile machine. 

As the load is increased, a stage is reached when the crack 

starts propagating. The rate of propagation is very slow and the 

strain output of the gauge is recorded with a strain indicator at 

regular intervals of time. Starting from zero when the load is 

zero, the strain increases as the crack comes near the strain 

gauge, reaching a peak. The strain then starts decreasing as the 

crack passes the strain gauge. Knowing the peak strain, the 

crack length, K is calculated making use of the numerical 

X o 

programme- 

For comparison purpbses, is determined by two other 

methods from the load recorded in the strip chart recorder 

corresponding to the peak strain. In the first method, comes 

out as a closed form solution described in section 2.B. 1 and in 

Appendix B and is called through approach' in this study. 

The second method finds K directly through the finite element 

1C 

scheme developed in chapter II. 
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Here it may be noted that the value obtained by finite 
element analysis is preferred to the one obtained through G 
approach as explained in Sec. 2.8.1. 

4.2 NUMERICAL SIMULATION 

The finite element programme developed in chapter II to 
determine strain field in a DCB specimen with thin cantilevers 
has been used for two purposes , (i) to develop a relation between 

and measured strain for a chosen crack length and (ii) to 
determine K in terms of load at cantilevers' ends. 

ICZ 

4.2.1 Relation between and Measured Strain 

A relationship is developed between the strain at a point 
where strain gauge is bonded and the corresponding The 

relationship is independent of the cantilever loads. 

It should be noted that the specimen is made of hardened 
alloy steel with very high yield stress. Consequently, the 
plastic zone at the crack tip is negligibly small especially for 
this study where loads and strain fields are small. The strain 
field can be considered to be elastic. Then for a given geometry 

and crack length strain c at any point of the specimen is 

proportional to the applied load, P. 

c a P ... (4.2) 

and k a P ...(4.3) 

By eliminating P, one obtains 

r oi e ... (4.4) 

'-J 

Thus any value of P can be chosen for running the computer 

programme. However, P should not be very large as it would 

introduce geometric nonlinearity. 
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For example with crack length a— 3£>.5 mm -for the specimen 
size given in chapter III, the programme was run -for arbitrarily 

chosen cantilever load P = 60 N. The peak strain £ along the 

P 

length of strain gauge is found to be related to K bv the 

I 

expressi on 


r: 0.00112 _ (4.5) 

Here K is expressed in MPa*/m and £: is in microns- This 

p 

relationship is used to process data- However, the relationship 
depends on crack length and whenever crack length is changed the 
programme has to be run to establish the rel ati onship- 

4.2-2 Numerical Determination of K 

1C 

The linear relationship between the load applied at the end 
of the cantilevers and K is obtained from the numerical 

I 

programme for a given crack length. Corresponding to the 
position of the peak strain, the crack length is known 
(Sec. 2.B.4). 


Considering the same example as in Bee. 4.2.1, the 
relationship is obtained as 

K 0-0395 P ... (4.6) 

1C 

Here also, the relationship is worked out by running the 
finite element code whenever crack length is changed. 

4.3 EXPERIMENTAL RESULTS 

A total of six experiments were conducted in this studv^. 
The crack length for each experiment is shown in column I of 
Table II- In each experiment, Load vs. Crack Opening 

Displacement (COD) and Strain vs. COD are recorded. The COD is 
directly recorded by the strip chart recorder. The chart speed 
and crosshead speed have been maintained at 10 mm/minute and 
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0. Imm/mi nute respectively. Hence 10 mm movement of the chart 
corresponds to a COD of 0.1 mm. The chart recorder plots the COD 
against the load applied. The strains measured through a strain 
indicator are then plotted on the same graph. 

For explaining the results of the experiment, the 

observations of the first experiment are considered. Figure 4.1 

shows two curves for experimental results. It is found that the 

observed peak strain is 138 /j. Using this value of e and 

Eq.4.4, the value of K is obtained as 1.55 MPaVm. This is 

referred to as experimental K in column IV of Table II. 

ic 

The load corresponding to the peak strain can also be read 
from the same figure as 5S N. Also, in this particular 
experiment, the strain gauge was at a distance of 38 mm from the 
hinges<load points). As mentioned in Sec. 2-8.4, the crack 
length corresponding to the peak is 36.5 mm. Now, for these 
values, the value is predicted as 2.29 Mpaym through 

approach using Eq. 2. 18. 


Similarly invoking linearity, K can also be calculated as 

xo 

1.89 HPaym using Eq. 4.6. The numerical value is less than the 
result obtained through Eq. 4-1 by 18 X. 


The difference in the numerical and experimental values is 
IB X. It is also seen that the value obtained from the 
experiment is smaller. 


Figures 4-2 through 4-6 give the graphs for the other five 

experiments conducted. In each experiment, measured strain 

increases as the crack tip approaches the strain gauge, reaching 

a sharp and distinct peak- In these experiments also, the stress 

intensity factors are calculated as explained above- The 

corresponding analytical and numerical values of also are 

calculated. The results are tabulated in Table II. The last 

column shows the percentage difference between experimentax anu 

numerical F values, 
xc 
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Table II shows that results of numerical simulation and G 

I 

approach match well. In all the six cases, G^ approach predicts 

K higher but within 18 X. This difference is expected as 
rc 

explained in Sec. 2.8.1 and so this validates the finite element 
programme. However the results of numerical simulation is 
regarded to be more accurate as explained in Sec. 2.8.1. 
Therefore the experimentally recorded is compared with 

obtained through numerical simulation in Col. VII of Table II. 


Measured strain is consistently lower. Accuracy of small 
strain gauges of 0.2 mm gauge length has been checked separately 
through a static test (Sec 3.6) and the error was found to be 
within 5 7.. Imperfect bonding with entrapped air bubbles is 
thought to be the reason for the discrepancy since such a 
situation will lead to change in strain field leading to 
incorrect analysis. 


It is expected that the strain gauge reading would be zero 
after the crack has passed the gauge since the location of the 
strain gauge is such that the bending strains are zero. Instead 
it is seen that the strain reading assumes another level and 
remains there even after unloading. The difference vanes 
considerably from one experiment to another. One of the specimen 
was checked after the two cantilevers were bonded under high 
pressure and temperature but before it was loaded. A slight 
cur-vature was observed which developed during the curing. As the 
strain gauge was glued to the prestressed DCB assembly, it did 
not show any strain until the crack passes through the strain 
gauge and prestressing is relieved. 


In order to keep the prestressing to low values, a new 
fixture was made after the first four experiments. The improved 
fixture provides ground flat and ver-> stiff surfaces above and 
below the specimen while curing. Experiments 5^ and 6 we^-e 

conducted on specimens prepared this wa> . Of the-e 



experiments, the effect is evident in the latter where the drop 
is only 10 microns compared to a peak of 120 microns. 

To sum up, a small strain gauge of 0.2 mm gauge length 

bonded at 45*^ to the crack plane is able to measure strain at 

about 2 mm radial distance from the crack tip. A distinct peak 

in measured strain is recorded as the crack passes under the 

strain gauge. Knowing the» peak strain, K is determined using 

xo 

the computer programme. For comparison, the load at the ends of 
the cantilever recorded at the instant of peak strain provides 
the through numerical simulation. 
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Fig. 4.6 Experimentally Recorded Strain and Load for 
Expt. No. 5, 


micro strain 
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CHAPTER V 

CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 


the work has been to experimentally determine the 
stress intensity factor by measuring strain near crack tip 
through a strain gauge in a double cantilever specimen with thin 
cantilevers. A finite element software has been developed for 
estimating the strain field which has been used in selecting a 
proper position and orientation for fixing a small strain gauge 
of 0.2mm. 


The finite element code was validated by solving two 
standard problems. The software was then used to determine the 
stress intensity factor as well as the strain field for the DCB 
specimen having thin cantilevers. The numerical value of stress 
intensity factor is in good agreement with the one obtained 
through energy approach. 

The DCB specimen has been made by bonding two thin strips of 
hardened alloy steel with epoxy. The curing was done at a 
temperature of 130*^ Celsius under a pressure of 1.2 MPa. While 
bonding, a precrack was introduced by keeping a BOPP sheet 
between the strips up to the required length. The strain gauge 
of gauge length 0.2 mm is bonded to the side face of the specimen 
at 45° with the help of a fixture made of self sticking paper. 

The specimen is loaded in mode I in an INSTRON machine. 
The experiment is conducted quasi -stati cal 1 y under controlled 
displacement conditions. As the crack grows and approaches the 
strain gauge, the output from the strain gauge increases till it 
reaches a peak. The peak corresponds to the location of the 
crack tip at about 1.5 mm prior to the strain gauge measured 
along the crack plane. Further crack extension reduces the 
strain. The measured peak strain gives the value of critical 
stress intensity factor through the relationship between strain 
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and SIF which is obtained through the computer programme. This 
experimentally obtained value o-f has been compared with the 
numerical prediction calculated using the load corresponding to 
the inst&nt at which strain reaches the peak value, with the help 
of the numerical programme. 

The measured strains are found to be lower by B to 55 % than 
the values predicted by the finite element analysis. The reason 
is thought to be imperfect bonding between the 
cantilevers due to entrapped air bubbles. This affects the 
stress/strai n field and make the analysis imperfect. 

In future, bonding of the specimen may be done in a vacuum 
bag which will avoid air entrapment and give a better bond. The 
method can be extended to dynamic problems and composite 
mater i al s. 
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APPENDIX A 

The behaviour of quarter point element is explained below 
£103. The derivation is for 1-D and can be extended to higher 
dimensions since the nodes on any sides of a rectangular element 
are not affected by those on the other three sides. 


An element having side nodes can be made to display a square 
root singularity in stress or strain by suitably defining its 
geometry. Consider for example the three noded bar element as 
shown in Fig. A. 1 (a). 

The displacement u at any point is given as 

u = Nu+Nu+Nu ... CA- 1) 

1 1 2 2 3 3 

and the coordinate x at any point is 

X = N X + N X + N X ... (A. 2) 

1 1 2 2 3 3 


12 3 

o o o 

2 f = -1 ?= 0 ?= +1 

Figure A. 1 (a) . 


1 

o- 


?= -1 ?= 0 
L^rL-.4-.- 
* — *■ X ,u ^ 


■(l-r)L- 


K= +1 


Figure A. 1 (b) . 


where N s are the shape functions defined as 

t 


N = - ? (1-^) 

12^ 

N = (1-?^) 

2 


N = 

3 


1 

2 


( 1 +^) 


. . . CA« 3) 
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Substituting in Eq.A.2, 


2 1 22 3 


. . . (A-4) 


If the third node is shifted to a point which is at a distance rL 
from node 1, it can be shown that r ^ - for the stress to be 
singular at node 1 and that the nature of singularity is a square 
root one (Fig A«l.(b))- From the figure, 

X = 0, X = rL- X == L 

12 ’3 

Hence from equation A- 4, 


= 5 [ 


r-2r?^ + ? + ?^ 


(A.5> 


Now dif ferentiating displacement u with respect to x , 


^ = in ^ 

dx dK 


[If] 


. . . (A. 6) 


where CJ3 denotes the Jacobian given by 


4 r ? + 1 + 2 ? 


(A. 7) 


Substituting in Eq.A.6, 


in = in 

Sk 


( -4 r ? + 


1 * 2 f > ] ‘ 


. . . (A. 8) 


At node 1, ? = -1 and K = 0- Therefore, 


Jli _ ^ 

dx ?=-i &K 


[k<^r + l- 2)j - ^ 


T< 4 r - 1 ? 
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For 


da 

dx t,=-± 


-+ ou 


-!i;(4r~l)=0 ^ r - - 

4 


(A. 9) 


atJEtituting in Eq.A-5, 




or , ( 1 4 ^ ) 


= J¥ 


. . . CA. 10) 


Substituting it in Eq.A-.7 along with the value of r, J simplifies 
to 

J _ ~ 1 4- ^ ) 


i . e. , J ’ = l]"^ ■ ■ • 

Using equations A. 3, A. 10 and A. 11 in A. 8, one obtains 


da 

dx 


[( ? ■ IcJ "C L ■ £ \x )“2 "C £■ - ££ J b 


which means stress is inversely proportional to square root of 
the distance' from the node. 
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APPENDIX B 

Equation used in calculating through has been cierived 
here. B can be calculated by the -formula CSD, 

P _ 

“21^ ...(B.l) 

i^'hore Ei is the thickness, a is the crack length and C is the 
compliance o-f the specimen and P is the applied load. C is 
de-fined as C = ^ where u is the displacement o-f the point loads. 

For a cantilever with length L, moment o-f inertia I and 

3 

Young's modulus E under the action of a -force P, u = In the 

case o-f a DCB specimen, the length of each of the cantilevers is 
equal to the crack length and since there are two cantilevers. 


C 




. . . (B.2) 


Elubsti tuting Eq- B.2 in B.l, the expression for B is obtained as 


G 

1 


12 P^a^ 
2 B a^ 




) 


3 Z 

Using the relation = — g. — (l~i^ > -for plane strain ^ the equation 
yi el ds 


i: 

1 


2 /3 P a 

.3 yZ . ^ 2 . 

B h (l“v ) 


. . . (B.4) 



66 


REFERENCES 


1. N. Takeda, R. L. Sierakowski, C. A. Ross Se L. E. Malvern 

-Dei ami nati on Crack Propagation in Bal i i st i cal 1 y 

Impacted Glass/Epoxy Composite Laminates. 

-Experimental Meehan ics ^ 22 <19S2^ pp 20-5 

2. Prashant Kumar, D. K. Sarkar & S. C. Gupta 

-Rolling Resistance of Elastic Wheels on Flat Surfaces 
-wea;-, 126 (1988} , pp 117-129 

3. A. J. Rosakis, J. Duffy & L. B. Freund 

-The Determination of Dynamic Fracture Toughness of AISI 
■^340 steel by the shadow spot method- 

— Journal of Phys, solidSf Vol. 32, Ho, 4 (1984), pp 443—6-0 

4. A- Shukla, B. D. Agarwal £c Bharat Bhushan 

-Determination of Stress Intensity Factor in Orthotropic 
Composite Materials using strain gauges. 

—Engineering Fracture Mechanics, Vol. 32, No, 3 (1989), 

pp 469-77 

5. J. R. Berger & J. W. Dally 

—A Spaciaily Overdetermi ned Analysis for Propagation 
Poughness Using Strain Gauges. 

—Mechanical Research Communications , Vol 17 (2), pp 9 s—r 

6. T. Nishioka, T. Murakami, H. Uchiyama, K. Sakakura & H. 
Kittaka 

—Specimen Size Effects on Dynamic Crack Propagation and 
Arrest in DCB Specimens. 

—Engineer ing Fracture Mechanics, Vol 39, No, 4 (1991), 

pp 757-767 

7. 0. Kolenik 

-On the Physical Meaning of the J-Aa 
-Engineer ing Fracture Mechanics, Foi 
pp 403-412 

8. David Broek 

-•Elementary Engineering Fracture Mechanics. 

—Martinus Nijhoff Publishers, 


Curves 

38, No, t (1991), 


9. R. D. Cook ^ 

-Concepts and Applications of Finite Element Anaiyexs. 

-John Hi ley & Sons, 


10 . 


R. S- Barsowm 

-On the Use of Isoparametric Finite Elements 
Fracture Mechanics. 

—International Journal of Numerical Methods in 
Vol, 10 (1970), pp 25-3-^ 


in Linear 


Engineering, 



67 


11.. Lm P. Harrop 

-The Optimum Size of Quarter — Point Crack^Tip Element. 
-International Journal of Numerical liethods in Engineer ing ^ 
Vol. 17, pp 1101-3 

12. Dr. K. Ramesh, Department of Hechanical 
Engineering, I I T, Kanpur. 

13. S. Q. zhang, B. Z. Jang, B. T. Valaire & J. C. Suhling 

-A New Criterion for Composite Material Mi.'jed Mode Fr a,?— t Lir e 
Analysis. 

-Engg. Fracture MechariicSf Vol 34, no, 3 (1989 pp 749-769 

14. E. P. Popov 
-Mechanics of Materials- 

-2nd Edition, Prentice Hall , Inc. 

15. M. D. Narayanan 

-M. Tech Thesis <1988). 



